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ABSTRACT 

We use hydrodynamical simulations from the OWLS project to investigate the dependence 
of the physical properties of galaxy populations at redshift 2 on the assumed star formation 
law, the equation of state imposed on the unresolved interstellar medium, the stellar initial 
mass function, the reionization history, and the assumed cosmology. This work complements 
that of Paper I, where we studied the effects of varying models for galactic winds driven by 
star formation and AGN. The normalisation of the matter power spectrum strongly affects the 
galaxy mass function, but has a relatively small effect on the physical properties of galaxies 
residing in haloes of a fixed mass. Reionization suppresses the stellar masses and gas fractions 
of low-mass galaxies, but by z — 2 the results are insensitive to the timing of reionization. The 
stellar initial mass function mainly determines the physical properties of galaxies through its 
effect on the efficiency of the feedback, while changes in the recycled mass and metal fractions 
play a smaller role. If we use a recipe for star formation that reproduces the observed star 
formation law independently of the assumed equation of state of the unresolved ISM, then 
the latter is unimportant. The star formation law, i.e. the gas consumption time scale as a 
function of surface density, determines the mass of dense, star-forming gas in galaxies, but 
affects neither the star formation rate nor the stellar mass. This can be understood in terms of 
self-regulation: the gas fraction adjusts until the outflow rate balances the inflow rate. 

Key words: cosmology: theory - galaxies: formation - galaxies: evolution - galaxies: fun- 
damental parameters - methods: numerical 



1 INTRODUCTION 

Understanding the star formation (SF) history of the Universe rep- 
resents one of the most fundamental pieces of the galaxy formation 
puzzle. There exists significant uncertainty as to which physical 
processes control the rate of star-formation in galaxies of a given 
mass. At the coarsest level, the rate at which gas enters galaxies is 
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controlled, at high redshift, by the rate of growth of the dark mat- 
ter halo, while at low redshift it also depen ds sensitively on the 
rate at which gas can cool i nto galaxies (e.g.lWhite & Rees|[l978l: 
i Hernquist & Springe j|2003l;IChoi & Nagaminell201filSchave et al.l 
|2010| ; Ivan de Voort et al.ll201 lh . However, matching observational 
constraints on galaxy masses, ages and star formation rates (SFRs) 
is much more complex than this simple picture would suggest, and 
although many attempts have been made to model the formation 
of a galaxy population using complex hydrodynamical simulations 
and semi-analytic models, there exists no clear consensus on what 
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combination of physical processes are required to explain the dis- 
tribution of galaxies seen in the local Universe. 

Difficulties in the simulations begin where gas becomes dense 
enough to form stars. This interstellar medium (ISM) gas has 
an extremely complex structure. Amongst other processes, mag- 
netic fields, turbulence, cosmic rays and radiative transfer may all 
play some role in determ ining the rate at which stars form (e.g. 
iMcKee & Ostrikeife007t) . For this reason, most cosmological sim- 
ulations treat the gas in the ISM with very simple "sub-grid" pre- 
scriptions. Additionally, the limited numerical resolution of the 
simulations prevents a detai led modelling o f SF, and they need to 
rely on em pirical laws (e.g . Kennicutt 1998). Howe ver, recent ob- 
servations dKennicutt et alj|2007l ; iBigiel et al1l2008l) show that the 
star formation rate surface density is a function of the molecular 
hydrogen density, or the surface density of c old gas dSchavell2004l : 
iRrumholz et al.ll201ll;lGlover & Clarldl2012l) . 

Furthermore, the stellar initial mass function (IMF) assumed 
in the simulations is not predicted from first principles, but is im- 
posed. Observational determinations of the stellar IMF are very 
challenging outside the solar neighbourhood so the local IMF is 
usually assumed to be universal and is applied to all galaxies at 
all redshifts. In addition to uncertainties related to processes that 
take place inside galaxies, the ultra-violet (UV) background pro- 
duced by the first generation of galaxies, ionizes the Universe and 
bathes gas in ionizing radiation. This ionizing UV background can 
then strongly suppress the infall of gas i n to low-mass haloes (e.g . 
Ouinn et all 1 19961: lOkamoto et alj [20o3 ; iPawlik & Schavd [20091 ; 
HambricketalfcOllI) . 

It is therefore important to assess, independently, how each of 
these uncertain physical processes affects t he propert i es of galax- 
ies formed in cosmological simulations. In lHaas et alj J201 3t here- 
after Paper I) of this series we use cosmological hydrodynamical 
simulations from the Overwhelmingly Large Simulations (OWLS; 
ISchave et alj|20 loT) to investigate the effects of cooling and feed- 
back on the galaxy population at z = 2. In this companion paper we 
turn our attention to other physica l processes. This wo rk (together 
with Paper I) complements that of ISchave et alj d2010t) . where the 
cosmic SF histories predicted by the OWLS simulations were ana- 
lyzed. 

In Paper I we describe the behaviour of the simulations in de- 
tail, and in particular we focus there on the effects of radiative cool- 
ing and energetic feedback from star formation and AGN. One of 
the main conclusion from Paper I is that the SF is self-regulated by 
the interplay between gas cooling onto galaxies and the feedback 
that offsets this gas accretion. The cooling rate onto galaxies is, for 
a fixed halo mass and redshift, mainly determined by the radiative 
cooling rate, which itself depends on chemical composition. We 
found that for many integrated galaxy properties the uncertainties 
in models for cooling and feedback give rise to a very large spread 
in galaxy physical properties (see Fig. 2). In this paper we discuss 
the importance of cosmology, reionization, the prescription for the 
ISM, the SF law and the stellar initial mass function. We will show 
that uncertainties in these processes are of secondary importance 
for the total amount of stars formed and for the SFR of simulated 
galaxies. However, we will also show that some of the processes 
discussed in this paper are important for the amount of mass resid- 
ing in the star forming ISM of galaxies, which we will argue can be 
explained in terms of self-regulated star formation. 

The structure of this paper is follows: in Sec.[2]we briefly de- 
scribe the simulations used in this study and the numerical tech- 
niques we employ. In Sec.|3]we describe how galaxy properties de- 
pend upon the physics included in the simulation. In this paper we 



Table 1. Overview of the values of the cosmological parameters according 
to WMAP3 {OWLS reference) and WMAP1 (as used in the Millennium 
Simulation). Symbols have their usual meaning. 
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consider the effects of: cosmology (Sec l3.lt . reionization (Sec. l3.2t . 
the effective equation of state of the ISM (Sec. 13.31 1, the star forma- 
tion law (Sec. [33), and the stellar IMF (Sec.[33J. Finally, in Sec.g] 
we summarize our findings and conclude. 



2 NUMERICAL TECHNIQUES 

T he simulations comp rising the OWLS project are described fully 
in ISchave et alj feoioh . Here we briefly summarize the reference 
simulation, with a focus on the physical prescriptions relevant for 
this paper. This simulation will be referenced throughout this paper 
as 'REF'. 

All simulations are performed with an extende d version of the 
N-Body Tree/SPH code Gadget3 (last described in Springel 2005) 
in periodic boxes of 25 co-moving /r'Mpc and contain 512 3 dark 
matter and the same number of baryonic particles (which can be 
either collisionless 'stars' or collisional 'gas' particles). The parti- 
cle mass of the simulations we use here is 8.68 x 10 6 M Q for dark 
matter and 1.85 x 10 6 M G for baryons (initially, the baryonic par- 
ticle masses change in the course of the simulation due to mass 
transfer from star particles to gas particles). The gravitational soft- 
ening length is initially fixed in co-moving coordinates at 1/25 the 
mean inter-particle spacing. Below z = 2.91 the softening is fixed in 
proper units, at 0.5 /j~'kpc. We provide tests showing convergence 
of our results with respect to simulation box size and the particle 
number in the appendix of Paper I. 

The cosmology assumed in the reference simulation is sum- 
mariz ed in Table [p and is deduced from the WMAP 3 year re- 
sults dSpergel et al7[2007l) . The re sults are largely cons istent with 
the more recent WMAP7 results dKomatsu et al .1 [2oTTh . the most 
notable difference is in cr 8 , which is 2.3 a lower in WMAP3 than 
in WMAP7. The primordial helium mass fraction is set to 0.248. 

As the subgrid model variation is the main power of the OWLS 
suite, we will now describe the parameters and subgrid models used 
in the reference simulation, before varying them in later sections. 
Radiative cooling and heating are calculated element-by-element 
by explicitly following the 1 1 elements H, He, C, N, O, Ne, Mg, 
Si, S, Ca and F e in the presence of the Cosmic Microwave Back- 
ground and the Haar dt & Madaul d200lh model for the UV/X-ray 
background ra d iation from quasars and galaxies, as described in 
IWiersma et all [2009a). The timed relea se of these ele ments by 
stars is followed as described by IWiersma et alj {2009b ) . The gas 
is assumed to be optically thin and in photo-ionization equilibrium. 
The simulations model hydrogen reionization by switching on the 
lHaardt & Madaul d200ll) background at z = 9. Helium reionization 
is modelled by heating the gas by an extra amount of 2 eV per 
atom. This heating takes place at z = 3.5, with the heating spread 
in redshift with a Gaussian filter with a(z) = 0.5. As shown by 
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Figure 1. A graphical representation of a galaxy in 10 12 ' 5 M© halo in 13 of our simulations at redshift 2. The colour coding shows the gas density in a slice 
of 100 h kpc thickness, divided by the mean density of the universe. All frames are of size 100 co-moving kpc/h and are centered on the position of the 
galaxy in the 'REF' simulation. All frames have a thickness of 100 co-moving kpc/h. The images are oriented so that the galaxy is face on in the 'REF' 
simulation. Although all simulations form a disk, it is clear that some of the physics prescriptions considered in this paper can have a significant effect on 
galaxy morphology. 



IWiersma et aL (2009b), the reionization prescription used in these 
simulations roughly matches the temperature history of the inter- 
galactic medium (IGM) inferred from observations by Scha ve et al.l 

fcooch . 

In the centres of haloes the density and pressure are so high, 
that the gas is expected to be in multiple phases, with cold and 
dense molecular clouds embedded in a warmer, more tenuous gas. 
This multi-phase structure is not resolved by our simulations (and 
the simulations lack the physics to describe these phases), so we im- 
pose a polytropic effective equation of state for particles with den- 
sities nu > 10 cm 4 . These particles are also assumed to be star 
forming, as th is is the density required to form a molecular phase 
We set the pressure of these particles to P oc p 7 eff 5 
where y e g is the polytropic index and p is the physical proper mass 
density of the gas. In order to prevent spurious fragmentation due 
to a lack of numerical resolution we set y e ff = 4/3, as then the ratio 
of the Jeans length to the SPH kernel and the Jean s mass are in- 
dependent of density ( Schave &DallaVecchi3l2008h . The normal- 



ization of the polytropic equation of state is such that for atomic 
gas with primordial composition, the energy per unit mass corre- 
sponds to 10 4 K, namely (P/k = 1.08 x 10 3 K ctrr 3 for ra H = 10 _1 
cm 3 ). The implementat i on of SF is stochastic, as described in 
ISchave & Dalla Vecchial d2008h . with a pressure-dependent SFR, 
obtained from the assumption of local hydrostatic equ ilibrium and 
the observed Kennicutt-Schmidt law (Kennicutt 1998). 

Energetic feedback from star formation is implemented kineti- 
cally, giving on average n = 2 of the SPH neighbours of each newly 
formed star particle a 'kick' such that the total energy in the out- 
flow corresponds to roughly 40% of the energy available from su- 
pernovae (SNe) of type II (i ncluding Ib,c). H ere rj is the mass load- 
ing factor. For our assumed lchabrierl d2003h IMF this corresponds 
to a w ind velocity of 600 km s -1 . See iDalla Vecchia & S chave 
( 2008) for more details on the kinetic implementation of SN feed- 
back. 

Haloes are identified using a Friends-of-Friends algorithm, as 
described in Sec. 2.2 of Paper I and we only show results over mass 
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ranges where the simulation achieves numerical convergence (Ap- 
pendix A of Paper I). This corresponds to haloes that contain at 
least 100 star particles when considering properties as a function of 
stellar mass and a minimum of 2000 dark matter particles when we 
plot properties against total halo mass. 



3 ISOLATING THE EFFECTS OF THE INPUT PHYSICS 

In this section we discuss each of the variations to the input physics 
in turn. Table[2]summarizes the simulations used in this paper and 
indicates in which subsection they are discussed. Bold-face val- 
ues indicate departures from the reference model. In the subse- 
quent sections we will discuss the sensitivity of galaxy properties 
to cosmology (Section l3.lt . reionization (Section 13.2b . the poly- 
tropic equation of state of the ISM (Section [33), the assumed SF 
law (SectionEH and the IMF (Section [33}. 

A graphical representation of the gas density of a z = 2 galaxy 
formed in the different simulations (except 'MILL' , see below) is 
shown in Fig.[T] The galaxy resides in a halo of total mass ~ 10 12 5 
M Q . It was first identified in the 'REF' simulation, where its posi- 
tion (defined as the centre of mass of all particles within 10% of the 
virial radius) was determined. The line of sight is along the z-axis, 
which is almost perfectly aligned with the angular momentum vec- 
tor of the gas within 10% of the virial radius (cos(</>) = 0.994). For 
the other simulations the image is centered on the same position, 
illustrating the remarkable similarity in the positions and orienta- 
tions of the galaxies. The 'MILL' simulation was run with differ- 
ent cosmological parameters, resulting in a different distribution of 
galaxies over the volume. This model is therefore not plotted in 
Fig-IU The physics variations discussed in this paper do not lead to 
large differences in the total stellar content of galaxies. Neverthe- 
less, we see in Fig.[T]that there are significant differences in galaxy 
morphologies, although a disk forms in all cases. These differences 
will be discussed throughout the remainder of this section. 

Fig. [2] shows, as a function of total halo mass the nine differ- 
ent galaxy properties we consider in this paper: medians of stel- 
lar mass fraction (/ star = M stm /M to[ , panel A), SFR (panel B), 
baryon fraction (/baryon = M baryon /M tot , panel C), fraction of star- 
forming gas (/ism = Mi SM /M tot , panel D) and gas mass fraction 
(/gas, halo = (M gas , , ot ai - M 1SM )/M t0U panel E). Then, as a function of 
stellar mass: medians of the molecular gas mass in the ISM (panel 
F), specific SFR (sSFR = SFR/M,, panel G), the inverse of the gas 
consumption timescale (SFR/Mism, panel H) and galaxy number 
density (the galaxy stellar mass function, panel I). I n panels F, G 
and I we compare t o obse rvatio ns of, respectivel y, Genzel et al.l 
d20ld) . IPaddietalJ d2007l) and iMarchesini etal] J2009h . as de- 
scribed in detail in Paper I. The black line in all panels is the 
'REF' simulation. The blue lines represent all physics variations 
discussed below. The red lines show the simulations that are dis- 
cussed in Paper I (the effects of metal-line cooling, and feedback 
from star formation and black hole growth). It is clear that the sim- 
ulations presented in Paper I strongly influence the amount of stars 
formed, the star formation rate and the gas and baryon fractions 
in the halo. Due to self-regulation, the physics variations discussed 
below mainly influence gas consumption time scales and ISM mass 
fractions. 



3.1 Cosmology 

To investigate the dependence of the galaxy properties on cosmol- 
ogy, and to facilitate comparisons to earlier work, we investigate the 



effects of changi ng the cosmolog ical parameters from the WMAP 
3-year results dSpergel et al]|2007l) used in the 'REF' simulation 
to the so-called 'concordance' or WMAP year 1 cosmology that 
was use d in many previous s tudies including the Millennium Sim- 
ulation ( Springel et "al] |2005l) . We will refer to this set of cosmo- 
logical parameters as the 'Millennium cosmology' and denote the 
model assuming this cosmology 'MILL' . The main differences be- 
tween the WMAP3 and Millennium cosmologies are firstly cr 8 , 
which is 0.74 in WMAP3 and 0.9 in the Millennium cosmology 
and secondly, the universal baryon fraction Q. h , which is 0.0418 in 
WMAP3 and 0.045 in the Millennium cosmology. Other parameter 
values are summarized in Table Q] The WMAP3 cosmological pa- 
rameters ar e largely consistent w ith the most recent 7-year results 
of WMAP dKomatsu etal .1201 ll) . although in WMAP7 the value of 
erg is slightly higher (0.8). 

Increasing <x 8 and Q. h leads to the formation of many more 
stars than in the 'REF' simulation. In order to roughly match the 
global SFR in the 'MILL' simulation, it is therefore necessary to 
double the amount of SN energy injected into the ISM per unit 
mass. The 'MILL' simulation therefore uses a mass loading of 
77 = 4 for SN-driven winds, rather than the rj = 2 used in the 
reference model. To isolate the effect of cosmology, we therefore 
compare 'MILL' to a simulation that uses the WMAP 3 cosmology, 
but the same SN feedback model as 'MILL'. We term this sim- 
ulation 'WML4' . In Fig. [3] the effect of the cosmology can thus be 
isolated by comparing the blue dashed and red dotted curves, which 
correspond to the WMAP3 ('WML4') and 'Millennium' ('MILL') 
cosmological parameters (as indicated in Table[T]l, respectively. 

The main effect of cr 8 is to set the time-scale for structure for- 
mati on. With a hi gher cr 8 structures of a given mass form earlier 
(e.g. |PeeblesHl993l) . The concentratio n of a DM halo at a given 
mass is set by formation time. Indeed, iDuffv et al] d2008l) showed 
that dark matter halo concentrations are significantly lower in the 
WMAP5 cosmology (cr 8 = 0.77) than in the WMAP1 cosmology 
(tr 8 = 0.90). Galaxy SFRs at a given halo mass could thus be influ- 
enced by the value of cr 8 through its effect on the halo gravitational 
potential. Comparing 'MILL' to 'WML4' in Fig. [3] the total stellar 
mass (panel A), the SFR (panel B) and the ISM mass (panel D) of 
high-mass haloes are all slightly higher for the WMAP1 cosmology 
than for the WMAP3 cosmology, due to the higher central densities 
in the WMAP1 cosmolo gy at a given halo m ass. Comparing these 
results to the global SFR JSchave et al. 2010), we note that cosmol- 
ogy has a much larger effect on the global SFR than on individual 
objects. This reflects the effect that the cosmological parameters 
have on the halo mass function. A higher tr 8 leads to more haloes 
of a given mass, and the larger number of haloes results in a much 
higher global SFR density. 

The primary effect of cosmology on the galaxy stellar mass 
function (panel I) is through its effect on the number density 
of haloes a t a giv en mass, i.e. the halo mass function, see also 
ICrain et alj (2009). At low masses, the M tM - / 8tar (panel A) and 
the M, - Mism (panel F) relations differ in 'MILL' and 'WML4' by 
~ 0. 1 dex). The stellar mass function at the low-mass end changes 
by up to 0.2 dex, because haloes of a given mass are more common 
in the simulation with the higher values of erg. Panel I also demon- 
strates that a higher cr 8 also leads to the formation of more massive 
galaxies (M sl a,- goes up to ~ 10 12 M Q for the 'MILL' simulation). 
We do not see an exponential cut-off in the predicted stellar mass 
function. This is because (as demonstrated by the gas consumption 
time-scale in panel H and the star forming gas mass in panel F), gas 
is being consumed rapidly in massive haloes where SN feedback is 
ineffective (see Paper I for more discussion on this point). Also, 
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Figure 2. Median relations between halo properties in all the simulations described in this work (blue lines) and in Paper I (red lines). The reference model is 
shown as a black curve in each panel. In Sec. 3 we consider subsets of simulations in more detail. In the first five panels we show, as a function of total halo 



Mbaryon/Mtot, panel C), fraction of star-forming 



mass, medians of stellar mass fraction (/ stiU = M stM /M tot , panel A), SFR (panel B), baryon fraction (/bar 
gas (/ism = A^ISM/Mot: panel D) and gas mass fraction (/ gaSi halo = (^gas, total _ ^lSM)/A/tot: panel E). The last four panels show, as a function of stellar mass, 
medians of the molecular gas mass in the ISM (panel F), specific SFR (sSFR = SFR/M,, panel G), the inverse of the gas consumption timescale (SFR/M[sm, 
panel H, where the two remarkable outliers are simulations with a more efficient star formation law, see Section l3~4l and galaxy number density (the galaxy 
stellar mass function, panel I). As described in the text, we show medians in bins along the horizontal axes for all haloes that satisfy the convergence criteria 
that apply to that specific panel. The horizontal, dashed line in panel (C) shows the universal baryon fraction for our chosen cosmology, the data points in 
panel (F) show a sub-set of the compilation studied in Genzel et al. (2010), the dotted black line in panel (G) shows the stellar mass - sSFR relation from the 
GOODS field (Daddi et al. 2007) and the shaded yellow region in panel (I) shows the galaxy stellar mass function of Marchesini et al. (2009). 
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Table 2. Overview of the simulations and the input physics that is varied in this paper. Bold face values represent departures from the reference model ('REF'). 
The first column gives the name of the simulation; the second column indicates the assumed cosmological parameters; the third lists the redshift of reionization; 
the fourth column specifies whether or not extra heat was injected during He reionization happens; the fifth column specifies the adiabatic index of the equation 
of sta te imposed on st ar forming gas; the sixth, seventh and eighth columns indicate the amplitude, slope and threshold of the SF law (K98 indicates the value 
from lKennicutilll998l) ; the ninth column specifies the IMF; the tenth column indicates whether the simulation includes SN feedback and the final column 
specifies the section in which the simulation is discussed. 
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WMAP1 
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WML4 
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1.4 (K98) 
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WMAP3 
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REIONZ06 


WMAP3 
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K98 


1.4 (K98) 
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Chabrier 
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REIONZ12 


WMAP3 


12 
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K98 


1.4 (K98) 


K98 


Chabrier 
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NOHeHEAT 


WMAP3 
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Chabrier 
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WMAP3 
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K98 
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Chabrier 
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3.3 
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WMAP3 
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5/3 


K98 
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Chabrier 


yes 


3.3 


SFAMPLx3 


WMAP3 
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yes 


3.4 


SFSLOPElp75 


WMAP3 
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4/3 
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3.4 


SFTHRESZ 


WMAP3 
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yes 


4/3 
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1.4 (K98) 


Z-dep. 
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3.4 


IMFSALP 


WMAP3 


9 
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4/3 


K98 


1.4 (K98) 


K98 


Salpeter 


yes 


3.5 



'radio-mode' AGN feedback, which is not considered here, is often 
assumed to cause this drop. 



3.2 Reionization 

The reionization of the Universe by quasars and galaxies has a pro- 
found effect on the temperatu re of the IGM and m ay impact the 
galaxy formation process (e.g. lOkamoto et a l]l2008t). Reionization 
is modelled in the simulations by turning on a UV+X-r ay back- 
ground from galaxies and quasars (Haardt & Madau 2001). All gas 
in the simulations is assumed to be optically thin to this radiation, 
which is assumed to be uniform and isotropic. In the 'REF' sim- 
ul ation, the redshift of r eionization is set to z, = 9. As shown 
in lWiersma et al.l J2009al) . the main effect of the UV background 
is to heat low-density, cold gas to T ~ 10 4 K. This leads to gas 
in haloes with virial temperatures < 10 4 K (M lot < 10 8 M o ) be- 
ing evaporated from their host haloes, and so effectively imposes a 
minimum halo mass in which star-form a tion can take place (e.g . 
Ouinn et all 1 19961: lOkamoto et alj |20o3 ; IPawlik & Schavd 120091 ; 
Hamb ricketal.ll201ll) . This mass scale corresponds to only ~ 10 
particles in our simulations. Hence, our simulations do not probe 
the regime where photo-heating is expected to be most important 
and they may underestimate the effect for somewhat higher masses 
because the effect of photo-heating on their progenitors is not re- 
solved. 

To isolate the effects of reionization on the galaxy population, 
we compare the 'REF' simulation to a simulation without reion- 
ization ('NOREION') and to two models in which the redshift of 
reionization is changed to z, ■ = 6 or z r = 12 from its default value 
z r = 9, which are denoted 'REIONZ06' and'REIONZW respec- 
tively. As described in Sec. [2] we match the observed temperature 
evolution of the IGM by modeling helium reionization as an addi- 
tional heat input of 2 eV per atom around z = 3.5. The simulation 
'NOHeHEAT' neglects this extra heating. 

Inspection of Fig. [4] (all panels) reveals that by z = 2 all 
knowledge of the redshift of reionization has been washed ou t 
of the resolved galaxy population, see also ISchave et alj d2010h : 



IWiersma etall j201ll) . Comparing 'NOREION' to the rest of the 
simulations, we see that the effects of neglecting the heating due to 
reionization are to increase the stellar masses (panel A), star form- 
ing gas masses (panel D) and gas masses (panels C and E), and 
hence the SFRs (panels B and G) of low-mass objects (M tot < 10" 
M Q ). Gas that should have been removed from these objects by the 
extra energy it receives from the photo-ionizing background is al- 
lowed to remain. This effect is also visible in the stellar mass func- 
tion (panel I), where at low masses (M star < 10 9 M Q ), a given stellar 
mass corresponds to a lower halo mass, making these objects more 
abundant. 

At high halo masses (M tot > 10" M Q ), neglecting reionization 
results in slightly lower halo gas fractions because more of the gas 
is in the ISM and in stars. 

We can examine the effect of helium reionization by com- 
paring 'REF' (solid, black curve) to 'NOHeHEAT' (dashed, pur- 
ple curve ). From all p a nels w e can immediately see that, consis- 
tent with lSchave et alj ( |2010|) . the extra heat input due to helium 
reionization has no effect on the properties of resolved haloes and 
galaxies. All massive objects have already formed a very signifi- 
cant proportion of their mass by the time that helium reionization 
takes place (z = 3.5) and this heat input has a negligible effect at 
gas densities typical of haloes. The unimportance of helium reion- 
ization holds for all of the properties of galaxies we investigate in 
this work. The extra heat input to the IGM from helium reioniza- 
tion is therefore only important for controlling the temperature of 
the low-density gas of the IGM iTheuns et al1l2002l) and has no ef- 
fect on the properties of galaxies formed in the simulations, at least 
for M, > 10 8 M o . 



3.3 The polytropic equation of state for high-density gas 

Large volume cosmological simulations lack both the resolution 
and the physics to model the multiphase ISM. We therefore im- 
pose an effective equation of state (EoS) for all gas particles with 
densities higher than «h = 0.1 cm 3 . The EoS we use is poly- 
tropic, P oc p>eff and in the 'REF' simulation y eff = 4/3. This 
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Figure 3. Same as Fig. [2] but only for the 'REF' simulation (black, solid curves), the 'MILL' simulation (red, dotted curves) and the 'WML4' simulation 
(blue, dashed curves). We show medians in bins along the horizontal axes for all haloes that satisfy the convergence criteria that apply to that specific panel 
(Appendix A of Paper I). The horizontal, dashed line in panel (C) shows the universal baryon fraction for our chosen cosmology in black, and the one for the 
Millennium cosmology in red. The 'WML4' (blue, dashed line) run uses the same cosmology as the 'REF' simulation, but the same feedback as the 'MILL' 
model. Therefore, to isolate the effect of the cosmology, the red, dotted curve should be compared to the blue, dashed curved. 



value of y c ff ensures that both the Jeans mass and the ratio of the 
Jeans lengt h and the kernel of the S PH particles are independent 
of density (Schaye & DaUa Vecchi3 l2008). provided that they are 
resolved at «h =0.1 cnT 3 , preventing spurious numerical fragmen- 
tation due to a lack of resolution. In order to isolate the effect of the 
EoS, we compare 'REF' to two simulations in which the slope of 
the EoS changed to be either isothermal, with y e ff = 1 C EOSlpO'), 
or adiabatic, with y c ff = 5/3 ('EOSlp67'). 

Changing the slope of the EoS has a significant effect on the 



visual appearance of the galaxy disk. Comparing 'EOSlpO' and 
'EOSlp67' in Fig.Q] we can immediately see that using a steeper 
EoS pressurizes the gas more strongly, resulting in a smoother gas 
distribution. Despite markedly changing the visual appearance of 
the disks, it is clear from Fig. [5] that the galaxy stellar properties 
(panels A and I) as well as the total amount of baryons associ- 
ated with the galaxies (panels C, D, E and F) and the SFRs (panels 
B and G) are virtually unaffected by the structure of the ISM on 
small scales. The only noticeable (although very weak) difference 
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Figure 4. As Fig. [5] but showing only the subset of simulations in which the reionization history is varied. In the 'REF' simulation (black, solid curve), 
reionization occurs at z = 9. The red, dotted curve shows a simulation without reionization ('NOREION'), whereas the blue, dashed and the green, dot-dashed 
curves have reionization redshifts of z = 6 and z = 12 respectively ('REIONZ06' and 'REIONZ12'). The magenta dot-dot-dot-dashed curve shows a simulation 
in which the extra heat input due to helium reionization at z = 3.5 is neglected ('NOHeHEAT'). The predictions only differ for the simulation that neglects 
reionization altogether. 



between the simulations is that in the very most massive galaxies, 
steepening the EoS leads to slightly more efficient SF at a given 
halo mass. This is because, at a given density, pressures (and hence 
star-formation efficiencies) are higher for a steeper EoS. 

It is striking that the properties of the galaxy are almost en- 
tirely insensitive to changes in the polytropic index (at least in 
the range y c g = 1 - 5/3 ). This is a consequence of ou r use of 
the prescription for SF of Schave & Dalla Vecchja] d2008l) . In this 
prescription the observed Kennicutt-Schmidt surface density law 



is analytically converted into a pressure law under the assump- 
tion of vertical, local hydrostat ic equilibrium. As demonstrated by 
ISchave & Dalla Vecchial §008 ) , this enables the same recipe to re- 
produce the observed SF law regardless of the equation of state, 
and without tuning any parameters. Because the pressure profile is 
determined by the gravitational potential, the SF we predict is the 
same if the potential is the same, irrespective of the assumed equa- 
tion of state (which does, however, change the thickness of the gas 
disk). 
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Figure 5. As Fig. [3] but showing only the subset of simulations in which the polytropic equation of state (EoS) that is imposed on high-density (n\\ > 0. 1 
cm -3 ) gas is varied. In the 'REF' simulation (black, solid curve) we use a polytropic EoS with a power law index of y e g = 4/3. The red, dotted curve shows a 
simulation with an isothermal EoS, y e ff = 1 ('EOSlpO'). The blue, dashed curve shows a simulation with an equation of state with an adiabatic EoS, y e ff = 5/3 
('E0Slp67'). The predictions are insensitive to the assumed equation of state because our prescription for star formation depends on the pressure rather than 
the volume density. 



Note that, as far as we know, all other groups use Schmidt-type 
SF laws, i.e. volume density SF laws, rather than pressure laws. In 
that case the SFR is expected to depend on the assumed equation of 
state, even if the gravitational potential and the pressure profile re- 
main unchanged. Because the observed SF law that the simulations 
are calibrated to match is a surface density rather than a volume 
density law, the calibration would have to be repeated if the equa- 
tion of state is changed. This is, however, generally not done, which 



means that the simulations can no longer be expected to reproduce 
the observed SF law when y c ff is changed. 

3.4 The star formation law 

Because large-volume cosmological simulations have neither the 
resolution nor the physics to model SF, the standard approach is to 
implement a simple, stochastic volume density SF law and to cal- 
ibrate it to match the observed Kennicutt-Schmidt surface density 
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Figure 6. As Fig. [5] but showing only the subset of simulations in which the parameters of the KS star formation law are changed. The 'REF' simulation 
(black, solid curve) uses a SF law that reproduces the observed KS law (Kennicutt 1998). The red, dotted curve shows a simulation where the amplitude of 
the KS law has been multiplied by a factor 3 ('SFAMPLx3'), and the blue, dashed line shows the effect of changing the slope of the KS law from 1.4 to 1.75 
('SFSLOPEIp75 7 ). These two models both have more efficient S F than 'REF' . The green, dot-dashed curve shows a simulation in which the SF threshold is 
a function of the gas metallicity ('SFTHRESHZ'), as predicted bv lSchavd <2004l) . The assumed star formation law affects the gas consumption time scale and 
the amount of star forming gas (ISM). It does not affect stellar masses or star formation rates, which can be understood if star formation is self-regulated (see 
text). 



law, 

1, =A(£ g /lM G pc- 2 )", (1) 

where n = 1.4 and A = 1.151 x 10 -4 M yr~' kpcr 2 a nd the law 
steep ens below a gas surface density X g ~ 10M Q pcT 2 jKennicuttl 
1998). In our simulations we do something similar, except that 
we implement SF as a pressure law, which has the advantage 



that no tuning is required, since we can convert the observed 
surface density law in a pressur e law under the assumption o f 
vertical hydrostatic equilibrium dSchave & Dalla Vecchial |2008). 
As we have seen in the previous section, another advantage of 
this approach is that the predicted SFRs become independent of 
the assumed equation of state of the star formin g gas (see also 
ISchave & Dalla Vecchi<Jl2008l ; ISchave et alj201fj|) . 
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In order to isolate the effect of changing the SF law, we com- 
pare the 'REF' simulation to three simulations in which the pa- 
rameters of the Kennicutt-Schmidt are varied. Simulation 'SFAM- 
PLx3' increases the normalization of the Kennicutt-Schmidt law, 
A, by a factor of three, which implies that for a gas particle at a 
given pressure, the SFR is a factor 3 higher in SFAMPLx3' than in 
'REF' . In simulation 'SFSLOPElp75' , the power-law slope of the 
KS-law is increased from n = 1.4 to 1.75. The normalisation of the 
Kennicutt-Schmidt law is chosen such that the SFR surface density 
is the same as 'REF' at S g = 1 M Q /pc 2 . As this is below the SF 
threshold, this KS-law forms stars more efficiently than the 'REF' 
simulation at all densities. The third variation on the 'REF' model, 
termed 'SFTHRESHZ' is a simulation in which the threshold den- 
sity for SF, n* n , depends on the gas metallicity, n* n oc Z -064 , such 
that the threshold density is equal to the reference simulation's SF 
threshold at a metallicity of 0. 1Z Q and there is a maximum allowed 
critical density of njjj = lOcirr 3 . This change mimics the metallic- 
ity dependence of the critical surfa ce density for t he formation of a 
cold, molecular phase predicted bv lSchavel 112004 1 . 

Panel (H) of Fig. [6] shows that, as expected, the gas 
consumption time scale is significantly (about a factor of 3) 
shorter for 'SFAMPLx3' and 'SFSLOPElp75' than for 'REF'. 
Model 'SFTHRESHZ' predicts nearly the same gas consump- 
tion time scale as 'REF' , although SF is slightly less efficient in 
'SFTHRESHZ' for all but the highest stellar masses. 

Intriguingly, the much shorter gas consumption time scale in 
models 'SFAMPLx3' and 'SFSLOPElp75' has virtually no effect 
on the galaxy stellar mas s (panels A and I) o r its SFR (panels B 
and G). This agrees with lSchave et all d2010l) . who found that all 
these simulations predict nearly the same cosmic SF histories. This, 
initially surprising, result can be understood by considering the dis- 
tribution of gas inside haloes. Panels (C) and (E) demonstrate that 
the total baryon fraction inside haloes does not depend on the as- 
sumed SF law, but that the amount of ISM gas, i.e. the amount of 
star forming gas, is very sensitive to the same parameters (panels D 
and F). Indeed, the response of the galaxies to a change in the gas 
consumption time scale, is to decrease the amount of star forming 
gas by a similar factor (com pare panels D and H ). This supports the 
explanation put forward by Scha ve et all d2010l) that star formation 
is self-regulated by the interplay between the available fuel supply 
and feedback processes. The feedback accompanying SF regulates 
the amount of gas condensing onto the ISM, an d hence SF itself, so 
that outflow roug hly balances inflow (see also Bouc he et al.ll201(il ; 
iDave et alj|2012ri . As the outflow rate depends on the SFR rather 
than on the gas consumption time scale, the SFR, and hence the 
stellar mass, remain the same if the SF law is changed. To produce 
the same SFR with a different SF law, the galaxy will have to adjust 
the amount of star forming gas. 

Naively one would expect that the metallicity-dependent SF 
threshold should have a large effect on the amount of stars formed 
in a simulation (zero metallicity gas has a SF threshold 100 times 
higher than in the 'REF' simulation), but Fig. [6] shows that this 
is not the case: the 'SFTHRESHZ' (green, dashed) curves closely 
track the 'REF' (black, solid) curved in all panels. This is another 
manifestation of the self-regulation described above. If the gas in a 
galaxy is insufficiently dense to form stars, then there are no pro- 
cesses to prevent gas from condensing into the galaxy and increas- 
ing the gas density until SF becomes possible. At this point the 
galaxy continues to form stars until the rate at which they inject 
energy into the galaxy is sufficient for the resulting galactic wind 
to b alance the rate of gas inflow . Seemingly contrary to our results, 
e.g. iGnedin & Kravtsovl d201 ll) found that a metallicity dependent 



cut-off is important for the SF properties of galaxies. The depen- 
dency on metallicity in those studies is, however, stronger than the 
one assumed here and the y focus on low-mass galaxies. 

Hopki ns"et al.l d201 ll) performed high-resolution galaxy simu- 
lations in which they varied the SF laws in a similar manner to that 
investigated in this section. Because of their much higher resolution 
(M& ~ 10 2 - 10 3 M ), which enables them to use a much higher 
SF threshold (« H = 100 cm -3 ), their feedback operates on a much 
smaller spatial scale. Because star particles form in the high-density 
peaks, their feedba ck is acting on much h igher-density gas than it is 
in our simulations. iHopkins et al.l a201 ll) find, in common with the 
results presented here, that the galaxy SFR is independent of the 
parameters of their SF law. However, in apparent contrast to our 
results, they find little dependence of the amount of gas in the ISM 
on these same parameters. This apparent di screpancy can be un- 
derstood by noting that in the simulations of Hopkins et al.l d201 lh 
there is a strong dependence of the high-density tail of the gas den- 
sity distribution function on the SF law parameters. They find that 
there is more gas at very high densities when the SF efficiency is 
lowered. The difference, therefore, results from the scale at which 
the feedback is imposed on the gas. For their simulations this occurs 
at much higher densities, and the feedback therefore only regulates 
the high-density end of the gas distribution - the gas where stars 
are actually forming. Our lower resolution simulations impose the 
feedback on larger scales and at densities comparable to our SF 
threshold, which coincides with our definition of the ISM (except 
for 'SFTHRESHZ'). Therefore, in both cases feedback from SF acts 
to regulate the amount of star-forming gas. As we do not a priori 
know on what scales feedback regulates SF in real galaxies, it is 
not clear which gas densities are sensitive to the SF law. 

3.5 The stellar initial mass function 

The assumed stellar IMF can potentially have a large effect on the 
galaxy population. The IMF determines the ratio of low- to high- 
mass stars. This, in turn, changes the integrated colours of stellar 
populations, their chemical yields and the number of SNe per unit 

stellar mass. 

In this section we compare 'REF', whic h assume s a Chabrierl 

(2003) IMF to a simulation that assumes a lSalpeterl dl955l) IMF 
('IMFSALP'). Above 1 M Q the two IMFs have similar shapes, but 
whereas the Salpeter IMF is a pure power law, the Chabrier IMF 
includes a log-normal decrease at the low-mass end, resulting in 
much lower mass-to-light ratios. Both IMFs are integrated over the 
same mass range (0.1 - 100 M G ). We increase the normalization 
of the SF law by a factor of 1.65 in the 'IMFSALP' simulation to 
take into account the difference in the number of ionizing photons 
produced per unit stellar mass (the observed SF law is based on the 
flux from massive stars combined with an assumed IMF). We do 
not, however, change the feedback parameters (wind velocity and 
mass loading are the same as in the reference model; for the effect 
of changing those parameters, see Paper I). Thus, in 'IMFSALP' , 
SNe inject 40% -1.65 ~ 66% of the total available SN energy into 
their surroundings (i.e. ~ 34% is assumed to be lost radiatively). 

Fig.[7]demonstrates that the total stellar masses of the galaxies 
(panel A), as well as the total amount of baryons in the halo (panels 
C and E), are not strongly affected by the choice of IMF. Inspection 
of Fig.Q]shows that the extent and morphology of the gaseous disk 
is also not strongly affected by the IMF. 

However, at a given halo or stellar mass, SFRs (panel B) are 
somewhat lower in the 'IMFSALP' simulation. This is because the 
Salpeter IMF produces less metals and returns less mass to the ISM 



© 201 1 RAS, MNRAS 000. [11041 



12 Marcel R. Haas et al. 





10.5 11.0 11.5 12.0 10.5 11.0 11.5 12.0 8 9 10 11 

Log 10 [M tol (M Q )] Log 10 [M tol (M )] Log 10 [M slQr (Mj] 




8 9 10 11 8 9 1011 891011 

Log 10 [M 5tar (M )] Log 10 [M star (M )] Log 10 [M slar (M )] 



Figure 7. As Fig. [3] but showing only the subset of simulations i n which the stell ar IMF is changed. The 'REF' simulation (black, solid curve) uses a[Chabriei 
1 2003) IMF, whereas 'IMFSALP' (red, dotted curve) assumes a lSalpeteil Jl955l) IMF. The star formation law has been rescaled for 'IMFSALP' to maintain 
consistency with the observations. This changes the gas consumption time scale (panel H) and, because of self-regulation, the ISM mass fraction (panels D 
and F). The other changes are small and can be attributed to the lower yields in 'IMFSALP' . 



than the Chabrier IMF. This decreases the available fuel supply 
both by keeping more mass locked up in stars and by decreasing 
the effect of metal-line cooling (see also Paper I). 

At a given mass, the amount of star-forming gas is signifi- 
cantly lower in 'IMFSALP' than in 'REF' (panels D and F). This 
occurs because the higher SF efficiency implied by the change in 
the SF law means that for a Salpeter IMF more stars are formed 
at a given gas density (see panel H). As discussed in Sec. l3.4l this 
means that the stellar population is injecting more energy into its 
surroundings at a given gas density (recall that the wind parame- 



ters are the same in the two simulations, only the SF efficiency has 
changed), and this allows the galaxy to self-regulate at lower gas 
densities. 



4 CONCLUSIONS 

We have analysed a large set of h igh-resolution cosm ological sim- 
ulations from the OWLS project ( Sch ave et al.ll201(jh . We focused 
on the baryonic properties of (friends-of-friends) haloes at redshift 
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2, while varying parameters in the sub-grid models for reionization, 
the pressure of the unresolved multiphase ISM, star formation (SF), 
the stellar initial mass function, as well as the cosmology. 

This paper complements Paper I, which focuses on variations 
in the energetic feedback from star formation and the effects of 
metal-line cooling and feedback from AGN. 

Our main conclusions can be summarised as follows: 

• The stellar mass, SFRs, and gas fractions of galaxies are in- 
sensitive to the stiffness of the equation of state that we impose on 
the un resolved, multiphase ISM. Thi s is a consequence of our use 
of the lSchave & Dalla Vecchk] |2008) recipe for SF, which ensures 
that the observed surface density law for SF is satisfied independent 
of the assumed equation of state. 

• The gas fraction of galaxies, i.e. the fraction of the baryonic 
mass in star-forming gas, is sensitive to the assumed SF law. If SF is 
assumed to be more efficient, i.e. if the gas consumption time scale 
at a fixed gas surface density is shorter, then the gas fraction will 
become lower. This is a result of self-regulation: the gas fraction 
increases until the formation rate of massive stars is sufficient to 
drive galactic winds that can balance the rate at which gas accretes 
onto the galaxies. As a consequence, the SFRs and stellar masses 
are insensitive to the assumed SF law. 

• In a cosmology with a higher erg structure formation happens 
earlier. Hence, galaxies residing in a halo of a fixed mass at a 
fixed time have somewhat higher stellar masses. The characteris- 
tic densities are also higher, which reflects the higher density of 
the Universe at the time the halo formed. These higher densities, in 
turn, cause feedback from SF to become inefficient at slightly lower 
masses if er 8 is higher. The differences in halo properties between 
different cosmologies are, however, much smaller th an the differ- 
ences between the cosmic SF histories we found in ISchave et al.l 
d2O10L see also [Springel & Hernquist 2003). This is because the 
halo mass function is sensitive to cosmology, which is more im- 
portant for the SF history than the relatively small change in the 
internal properties of the galaxies at a fixed time and halo mass. 

• For low-mass haloes (M < 10" M G , M, < 10 9 M Q ) the re- 
heating associated with reionization is important, although by z = 2 
the results are insensitive to the redshift at which reionization hap- 
pened, at least as long a s it happened n o later than z = 6, as 
required by observations jFan et all [2006). Without reionization, 
these haloes would host higher-mass galaxies with higher gas frac- 
tions. 

• The IMF mainly affects the physical properties of galaxies be- 
cause it determines the amount of energy and momentum that mas- 
sive stars can inject into the ISM per unit stellar mass formed. If this 
feedback efficiency is kept fixed when the IMF is changed, then the 
assumed SF law must change if it is to remain consistent with ob- 
servations. This will in turn modify the gas fractions, as explained 
above. The IMF also changes the metal yields and the rate at which 
mass is returned from stars to the ISM (recycling), but these effects 
play a smaller albeit still significant role. 

We conclude, also on the basis of the results presented in Paper 
I, that the integrated physical properties of galaxies are mainly de- 
termined by the efficiency of the feedback from star formation and 
AGN, while the efficiency of star formation (i.e. the star formation 
law) determines the gas fractions. To make further progress, it is 
therefore most important to improve our understanding of galactic 
winds. 
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